Dynamics of a driven probe molecule in a liquid monolayer 
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We study dynamics of a probe molecule, driven by an external constant force in a 
liquid monolayer on top of solid surface. In terms of a microscopic, mean-field-type 
approach, we calculate the terminal velocity of the probe molecule. This allows 
us to establish the analog of the Stokes formula, in which the friction coefficient is 
interpreted in terms of the microscopic parameters characterizing the system. We 
also determine the distribution of the monolayer particles as seen from the stationary 
moving probe molecule and estimate the self-diffusion coefficient for diffusion in a 
liquid monolayer. 
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Thin liquid films adsorbed on solid surfaces or confined in a narrow space between two 
solids represent a remarkable example of a two-phase physical system, in which an intrinsi- 
cally disordered liquid phase is spanned by and contends with the ordering potential of the 
solid. It is well appreciated by now that the behaviour of such films is markedly different of 
the customary behaviour of the bulk liquids: experiments reveal effects of solid-like or glassy 
response to an external shear, sharp increase of relaxation times and of the films' viscosity 
The departure from the bulk behaviour is progressively more pronounced the thinner 
the film is and is interpreted usually as the manifestation of intricate cooperative effects, 
whose nature is not yet completely elucidated (see, e.g. 0). 

In this Letter we study cooperative behaviour in a liquid monolayer on solid surface, 
arising in response to internal perturbances created by a single probe molecule (say, a charge 
carrier), which is driven by constant external force. We propose here a simple analytically 
soluble model, which allows us to determine the stationary, non-homogeneous distribution 
of the monolayer particles as seen from the probe molecule, and to derive a closed-form 
non-linear equation relating the terminal velocity of the probe to the magnitude of the 
driving force F. We find that in the limit of small F this equation reduces to the Stokes 
viscous-flow law, in which the friction coefficient £ is determined explicitly, in terms of the 
microscopic parameters characterizing the system under study. We show that £ combines 
the contributions of two different " dissipative" processes - temporary trapping of the probe 
by the potential wells, due to the interactions with the solid atoms, and the cooperative 
effects associated with the formation of a non-homogeneous distribution of the monolayer 
particles around the probe. The knowledge of £ also allows us to estimate the self-diffusion 
coefficient for diffusion in a liquid monolayer. 

We start with the formulation of our model. Consider a 2D homogeneous monolayer of 
identical interacting particles moving randomly on top of an ideal crystalline surface (Fig.l). 
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We suppose first that the particles' motion is induced by chaotic vibrations of the solid atoms. 
Second, we assume that the particle-particle interactions (PPI) are essentially weaker than 
the particle-solid interactions (PSI), such that the former don't perturb the wafer-like array 
of potential wells created due to the PSI (Fig.l). This permits us to adopt the conventional 
picture of the dynamics under such conditions (see and references therein): the particles 
migration proceeds by thermally activated rare events of hopping from one well to another 
in its neighborhood. The hopping events are separated by the time interval r*, which is the 
time a given particle typically spends in each well vibrating around its minimum; r* is related 
to the temperature and the barrier of the PSI by the Arrhenius formula. Further on, 
the PPI couple the dynamics of a given particle to motions of all others. That is, a particle 
escaping from a given well follows preferentially the gradient of the PPI potential (see, e.g. 
ITJII for discussion). The realistic PPI are repulsive at short scales and show a weak, long- 
range attraction at longer interparticle separations. Here we simplify the actual behaviour by 
neglecting, first, the particle-particle attractions and, second, by approximating the repulsive 
part of the interaction potential by an abrupt hard wall, imposed at the distance equal to 
the interwell spacing. In such an idealized model the hopping probabilities are naturally 
decoupled of the particles distribution and all hopping directions are equaly probable. The 
hard-core repulsion prevents multiple occupancy of any potential well - we suppose that a 
particle may hop into the target well only in case when the latter is unoccupied at this 
moment of time; otherwise, the particle attempting to hop is repelled back to its position. 

The key aspect in our model is the probe molecule (abbreviated as the PM in the fol- 
lowing), which is subject to an external constant (small) force F directed along the X-axis 
(Fig.l). We suppose that the PM also performs an activated hopping motion, constrained 
by the hard-core interactions with the monolayer particles. The typical time which the PM 
spends being trapped by each well is denoted as r, which is dependent on the probe-solid in- 
teractions and hence may differ from r*. Varying r we can mimic different possible regimes. 
For instance, the choice r = corresponds to the situation when the PM simply slides 
regardless of the surface corrugation. Now, in contrast to the monolayer particles the PM 
dynamics is anisotropic due to the applied force; that is, on escaping from the wells, the 
PM attempts to hop preferentially in the direction of F. Supposing for simplicity that the 
wells form a square (X, Y)-lattice of spacing a, we define the probability q(X,Y\X' ,Y') of 
an attempt to hop from the well at (X, Y) to the adjacent well at (X', Y') as: 

q(X,Y\X',Y') = Z- 1 exp((3F(X' - X)/2), Z = 4 cosh 2 (/3aF/4) (1) 

To take into account the hard-core repulsion between the PM and the monolayer particles 
we stipulate that the PM hops into the target well only in the case when it is vacant at this 
moment of time. Otherwise, it remains at its position. 

Let P(X,Y) (p(X,Y)) denote the probability that at time t the PM (a monolayer par- 
ticle) occupies the well with the coordinates (X,Y). To derive the closed-form equations 
describing the time evolution of these properties, we assume that correlations between the 
occupations of different wells decouple. It was shown recently in || that such a decoupling 
provides an adequate description of the driven particle dynamics in a one-dimensional sym- 
metric hard-core lattice gas, in which case the hindering effect of the gas particles is much 
more strong than in 2D. We thus expect that in the 2D situation under study such an ap- 
proach will yield correct dependence of the PM terminal velocity on the system parameters; 
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we do not expect, of course, to obtain exact numerical factors. Decoupling the correlations, 
we have 

4r* ^iXl = £ {_ (1 - p(X'X))p(X,Y) + (1 - p(X,Y))p(X'X)}, (2) 
at (X',y) 

which holds for all wells excluding the wells adjacent to the PM position. Next, we find 
r dP ^ Y) = £ {- q(X,Y\X',T) P{X,Y) (1 - p(X' ,Y')) + 

at (X',Y') 

+ q(X>, Y'\X, Y) (1 - p(X, Y)) P(X', Y')} (3) 

The sums in eqs.(0) and ([$]) run over all wells (X', Y') adjacent to the well at position (X, Y). 
Multiplying both sides of eq.(Q) by X and summing over all wells we obtain 

V = -(q + - q_) - - J2P(X,Y)(q + p(X + a,Y) - q_ p(X - a,Y)) = 
r T XY 



= ~ [q+ (1 - MAi = a, A 2 = 0)) - g_ (1 - / v (A a = -a, A 2 = 0))], (4) 

T 

where V is the X-component of the PM velocity, q± = exp(±/?aF/2)/Z and /v(Ai,A2) 
denotes the time-dependent probe-particle correlation function; the subscript "V" signifies 
that this property is itself dependent on the PM velocity. 

Dynamics of /y(Ai, A 2 ) is a superposition of two different processes - the particles diffu- 
sion and the motion of the PM. From eqs.(|2|) and (^) we find that in the continuous-space 
limit /y(Ai, A 2 ) obeys (see for details of derivation in the ID case): 



-MA 1; A 2 ) =£>*(_ + —)f v (\ 1} \ 2 ) + V — MA 1; A 2 ), ( - r ) ) 

where D* = o 1 /At* is the "bare" diffusion coefficient of an isolated particle. 

Let us discuss now the initial and boundary conditions to eq.(|5]). First, we suppose 

/v(Ai,A 2 )| AltA2 _ ±00 = MAi,A 2 )| t=0 = p, (6) 

which mean that correlations between P(X, Y; t) and p(X + Ai, Y + A 2 ) vanish in the limit 
Ai,A 2 — ► ±oo and that initially the particles were uniformly distributed in the potential 
wells with mean density p. Another set of four equations can be derived by analyzing the 
evolution of /v(Ai,A 2 ) in the wells adjacent to the PM position (see for more details 0). 
This gives the zero-flux boundary condition for Ai = 0, A 2 = ±cr and 



^ /v(Ai,0) 



<9Ai 



= -^M±M) (7) 

Al=±(T 
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Solution of the mixed boundary value problem in eqs. fl5|) to ([?]) yields the following result 
for the particle distribution around the stationary moving PM: 



f n M 9 , PA,, r cosh(P)-2S_ - 4^P) X l Pi 

A 2 ) = P + 2 p ex P (-— ) { S+ + 2S _ ^(-) KU-) ~ 

cosh(P)-2S - I 2n -i(P) T A Pr 
" 2 (1 + 4(S + + 2S_) } £ KUP) 12 - 1 ^ 1 n_ ~o~ (8) 

where S± are given by 

oo 77 / p\ co 

S + = - E T7^#7Wu and 3- = ^n-i(P) {ln[K 2n ^(P)]y (9) 



{Zn[# 2n (P)]}' 



71=1 



In eqs. (|) and @ P is the "Peclet number", P = VooCr/2D*, V'oo is the PM terminal velocity, 



the prime stands for the derivative with respect to P, r = y X 2 + A|, while K n (x), I n (x) and 
T n (z) denote the modified Bessel functions and the Tchebyscheff polynomials respectively. 
Inserting eqs.@ and @ into eq.(|) we arrive at a closed with respect to Voo equation. 
This non-linear transcendental equation is not analytically soluble, and we will thus seek an 
approximate solution, assuming that is sufficiently small, such that P < 1. Expanding 
/V(±<r, 0) to the second order in powers of P we obtain 



(g+ - - p) 

(r/a) + (q + + q-)ap/D* 



x [i + ^- ( ^7 2r )(1 : p) f ) 2 M ^" g : f )(1 : p) f ) + 7 + 1/2}], do) 

where 7 = 0.577... is the Euler constant. 

Now, several comments on the result in eq.(|T0|) are in order. We find, first, that P < 1 and 
our approximation is justified when the inequality — g_)(l— p)a 2 < 2D*r + 2(g + +g_)cx 2 p 
is fulfilled. This happens, namely, in the case when the applied force is sufficiently small, 
or when either p > 1/3, or D*r/a 2 > 1. We note also that in the limit P < 1 the second 
term in brackets in eq. (|It]) is much smaller than unity and can be safely neglected. Further 
on, eq.(|T0D shows that the PM velocity is controlled by the combination of two different 
factors - trapping by the wells and the hindering effect of the monolayer particles. The 
first factor dominates in the limit when D*r/a 2 ^> p, which behaviour may take place in 
the case when the barrier of the probe-solid interactions is much stronger than that for the 
PSI. In this limit eq.fllOD reduces to the essentially mean-field result ~ (q + — qJ)a/T u , 
where r u =t/(1— p). We remark that since r is appropriate for the motion in the stagnant 
layer, rather than in the bulk phase, we may expect that even in this limit of the well-stirred 
monolayer the terminal velocity will be much less than the corresponding (to this value of 
F) velocity in the bulk liquid. Next, the prevalence of the cooperative effects appears in the 
limit when the ratio D*t/<j 2 is less than p. Here we have ~ (q + —qJ)(l—p)D*/<jp, i.e. the 
terminal velocity is proportional to the "bare" diffusion constant of the monolayer particles 
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and vanishes when D* — > 0. This means that the particles have to diffuse away of the PM in 
order it can move. As a matter of fact, such a dependence mirrors an even stronger effect: 
the monolayer particles, whose motion is random with no preferential direction, accumulate 
in front of the obstacle (the PM), which moves at a constant velocity. 

Next, using the notations of eq.(|I|) and supposing that (3aF < 1 we cast eq.(10) into 
the following physically revealing form 

which can be thought of as the analog of the Stokes formula for the system under study; 
the friction coefficient £ is expressed through the microscopic parameters and displays the 
combined hindering effect of the potential wells and of the monolayer particles on the PM 
dynamics. The friction coefficient is proportional to r, which is related to the probe-solid 
interactions, and diverges when either D* — > or p — > 1. 

Further on, from the Einstein relation between the mobility and the diffusivity of a test 
particle (the PM) we find 

Dpm = (12) 
D* + 2 p D y ! 

where D = a 2 /At and Dpm are the PM "bare" and self-diffusion coefficients respectively. 

Let us now briefly discuss the density distribution of the monolayer particles as seen from 
the stationary moving PM. The local density in the nearest to the PM wells /y(±cr, 0) ~ 
p(l ± Vooa/D*), i.e. in front of (past) the PM the local density exceeds (is less than) the 
average density of the monolayer particles. Further on, at large separations of the PM, such 
that Ai A = D*/Voo, the density profile is characterized by 

MA 1; 0) « p (1 + (^) 2 (^V 2 exp(-A 1 /A)), (13) 

which means that the excess density in front of the PM vanishes exponentially with the 
distance. The parameter A can be identified as the characteristic length of the condensed 
region in front of the PM. We note, however, that this parameter is not very informative; 
particularly, it diverges when either D* — > oo or p — > 1, which is a bit misleading since in 
this limit the excess density tends to zero, i.e. the monolayer is homogeneous. A better 
parameter is the integral excess density Q along the Ai-axis, 

n = rd\ 1 (f v (\ 1 ,o) - p) « ~ g - )p( * ~ f f {-hp/2) + i- 7 } (i4) 

J a D*r + (q + + q-)a i p 



Eq.(14) shows that the condensed region is actually absent when either p— > 1 or /)* - >oo 
and increases when the driving force is increased. We also remark that when the parameter 
D*r/a 2 is sufficiently large, Q is a non-monotoneous function of p. 

Finally, we find that at large distances past the PM the density follows 

MA 1; o) « P (i + ^) 4 1/2 )' N» A ( 15 ) 
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Eq.(|15D displays two remarkable features: First, the particle density past the PM approaches 
its average value only as a power-law, which mirrors strong memory effects of the medium 
and signifies that diffusive homogenization of the monolayer, constrained by the hard-core 
interactions, is a very slow process. Second, it shows that the density is a non-monotoneous 
function of Ai; /y(Ai = —a, 0) < p and approaches p from above when Ai — > — oo. 

To summarize, we have studied dynamics of a probe molecule driven by an external 
constant force in a 2D monolayer of particles moving randomly on top of solid surface. We 
have proposed a microscopic model description of such a system evolution and calculated the 
PM terminal velocity V^,, self-diffusion coefficient Dpm and the stationary distribution of the 
monolayer particles as seen from the probe molecule. We have shown that this distribution is 
strongly asymmetric: Past the PM, the local density is lower than the average density p and 
tends to p as a power-law of the distance, revealing strong memory effects of the medium. In 
front of the PM the density is higher than the average - the monolayer particles accumulate 
in front of the PM, creating a sort of a "traffic jam" which impedes its motion. We show 
that the PM terminal velocity is controlled by the size of the jammed region, which is, in 
turn, dependent on the driving force, as well as on the density of the monolayer and on the 
rate at which the monolayer particles may diffuse away of the PM. The balance between the 
driving force and the rate of the monolayer homogenization manifests itself as a medium- 
induced frictional force exerted on the PM, which shows for small driving forces a viscous-like 
behaviour. This resembles in a way the situation described in |1(J, which work considered 
the dynamics of a single charged particle moving at a constant velocity above the surface 
of a dielectric fluid. Here the cooperative behaviour emerges because of the electrostatic 
interactions between the particle and the fluid molecules. The particle perturbs the fluid 
surface creating a "bump", which travels together with the particle increasing effectively its 
mass and thus producing a frictional drag force exerted on the particle. 
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cussions. Financial support from the COST Project D5/0003/95 and the FNRS is gratefully 
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Fig.l. A schematic picture of a liquid monolayer on solid surface. The black circle denotes 
the probe molecule and smaller grey circles represent the monolayer particles. Wavy lines 
depict the wafer-like potential landscape created by the particle-solid interactions. 
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